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ABSTRACT 

We present rotational velocities for individual components of eleven very low mass (VLM) binaries 
with spectral types between M7 and L7.5. These results are based on observations taken with the 
near-infrared spectrograph, NIRSPEC, and the Keck II laser guide star adaptive optics (LGS AO) 
system. We find that the observed sources tend to be rapid rotators (usini > 10 km s _1 ), consistent 
with previous seeing-limited measurements of VLM objects. The two sources with the largest wsini, 
LP 349-25B and HD 130948C, are rotating at ^30% of their break up speed, and are among the most 
rapidly rotating VLM objects known. Furthermore, five binary systems, all with orbital semi-major 
axes <3.5 AU, have component usini values that differ by greater than 3<r. To bring the binary 
components with discrepant rotational velocities into agreement would require the rotational axes to 
be inclined with respect to each other, and that at least one component is inclined with respect to the 
orbital plane. Alternatively, each component could be rotating at a different rate, even though they 
have similar spectral types. Both differing rotational velocities and inclinations have implications for 
binary star formation and evolution. We also investigate possible dynamical evolution in the triple 
system HD 130948A-BC. The close binary brown dwarfs B and C have significantly different vsini 
values. We demonstrate that components B and C could have been torqued into misalignment by the 
primary star, A, via orbital precession. Such a scenario can also be applied to another triple system 
in our sample, GJ 569A-Bab. Interactions such as these may play an important role in the dynamical 
evolution of very low mass binaries. Finally, we note that two of the binaries with large differences in 
component usini, LP 349-25AB and 2MASS 0746+20AB, are also known radio sources. 
Keywords: stars: binaries, visual; stars: low-mass, brown dwarfs; stars: fundamental parameters; 
stars: rotation 



1. INTRODUCTION 

Rotational velocity is an important diagnostic parame- 
ter for stellar objects, offering a window into the angular 
momentum evolution of a given source. A star's rota- 
tion can provide important clues to its formation and 
can furnish diagnostics of its interior structure and evo- 
lution. For instance, measurements of rotational velocity 
have been shown to correlate strongly with stellar activ- 
ity, possibly driving the magnetic dynamo responsible 
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for generating this activity ([Browning! [2008) . In addi- 
tion, rotational velocities have been shown to correlate 
with the ag e of a system, offerin g a tool for estimating 
stellar ages ([Delfosse et al.lll99l . 

The rotational behavior of very low mass (VLM) 
stars and brown dwarfs has been st udied by a num- 
ber o f authors in rec e nt years (Mohantv &: Basril 
20031 iBailer-JonesI [20041: IZapatero Osorio et al.1 120061: 
Reiners fc Basril 120081 120101 iBlake et al.1 120101 ) It has 
been shown that the brown dwarfs tend to be rapid ro- 
tators, and that the mini mum rotation rate is a func- 
tion of spectral type (i.e., IZapatero Osorio et al.l 120061 : 
iReiners fc Basril I2008T) . It has also been determined 
that the rotational velocities of brown dwarfs correlate 
with age, with VLM objects having very long spindown 
timescales, and that their rotational evolution is prob- 
ably dominated primarily by magnetic brakin g (e.g., 
IReiners fc Basril l200l IScholz et alJ[2009l How- 
ever, it appears that the activity-rotation relationship 
that is very strong among st M dwarfs tends to b reak 
down at these low masses (|Mohantv fc Basril [200l . In 
spite of this, activity in the form of radio emiss ion has 
been observed in a num ber of VLM systems ( e.g, iBergerl 
[20061 lOstenlt~aT1[2006l : lHallinan et aLll2008T ). In addi- 
tion, it has been shown that this drop in activity does not 
seem to be due to a reduction in magnetic field strengths 
in late M dwarfs, but might instead be due to their 
reduced temperature and hence reduced fractional ion- 
ization of their atmospher es (e.g.. IReiners fc Basril[2007l : 
lHallinan et aT1[2(M l2008h . 

The majority of previous studies have been performed 
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Table 1 

VLM Binary Sample 



Source Name 


RA 


Dec 


Estimated 


Discovery 


2MASS 




(J2000) 


(J2000) 


Sp Types a 


Reference 


K Band Mag. 


LP 349-25AB 


00 27 55.93 


+22 19 32.8 


M8+M9 


1 


9.569 ± 0.017 


LP 415-20AB 


04 21 49.0 


+19 29 10 


M7+M9.5 


2 


11.668 ± 0.020 


2MASS J07464256+2000321AB 


07 46 42.5 


+20 00 32 


L0+L1.5 


3 


10.468 ± 0.022 


GJ 569Bab 


14 54 29.0 


+16 06 05 


M8.5+M9 


■1 


-9.8 


LHS 2397aAB 


11 21 49.25 


-13 13 08.4 


M8+L7.5 


5 


10.735 ± 0.023 


2MASS J14263161+1557012AB 


14 26 31.62 


+15 57 01.3 


M8.5+L1 


6 


11.731 ± 0.018 


HD 130948BC 


14 50 15.81 


+23 54 42.6 


L4+L4 


7 


~11.0 


2MASS J17501291+4424043AB 


17 50 12.91 


+44 24 04.3 


M7.5+L0 


2 


11.768 ± 0.017 


2MASS J18470342+5522433AB 


18 47 03.42 


+55 22 43.3 


M7+M7.5 


8 


10.901 ± 0.020 


2MASS J21402931+1625183AB 


21 40 29.32 


+16 25 18.3 


M8.5+L2 


6 


11.826 ± 0.031 


2MASS J22062280-2047058AB 


22 06 22.80 


-20 47 05.9 


M8+M8 


6 


11.315 ± 0.027 



Note. — References - (1) Forveille et al. 2005 (2) Siegler et al. 2003 (3) Reid et al. 2001 (4) Martin ot al. 
2000 (5) Freed et al. 2003 (6) Close et al. 2003 (7) Potter et al. (2002) (8) Siegler et al. 2005 
a From discovery reference 



with seeing-limited observations, and most sources tar- 
geted are thought to be single. Known binaries have 
been included in various samples, and their rotational 
velocities have been derived from the combined light of 
both components. The rotational velocities of individ- 
ual binary components can potentially provide a unique 
look at the rotational evolution of VLM objects. If 
any differences are seen between the velocities of the 
binary components, it could have implications for the 
way in which these binaries formed, their early accre- 
tion history, or the operation of magnetic braking as 
a function of mass. For instance, iReiners et al.l (|2007f ) 
found that the components of the triple system LHS 
1070 (spectral types M5.5, M9, and M9) had differing 
wsini, with the higher mass component rotating about 
a factor of two more slowly than the two lower mass ob- 
jects. This allowed the authors to put constraints on 
the form of rotational braking in t he VL M regime. Fur- 
ther, iGomez Maaueo Chew et all (|2009[) found that the 
components of the young eclipsing binary brown dwarf 
2MASS0535-05AB (spectral types M6.5) have different 
rotational periods, with the primary component rotating 
more rapidly than the secondary. 

If the orbits of these binaries are known, rotational ve- 
locities provide a way to test the assumption that spin 
axes are generally perpendicular to the orbital plane. 
Work bv lHald (|1994[) found that, in general, binaries with 
separations <30-40 AU should have spin axes perpendic- 
ular to their orbital plane. However, some very close 
semi-major axis < -3 AU) binaries such as DI Herculis 
Albrecht et al|[2009D have been shown to have extremely 
misaligned axes. In contrast, m ore recent work on a very 
similar binary system (NY Cep. lAlbrecht et al.ll201lh has 
revealed no such misalignment, implying that the cause 
of the inclined axes is non-universal. It is also important 
to explore systems with wider separations that may not 
be subject to the same extreme tidal interactions as very 
close binaries. This has been done in the case of some 
T Tauri binary systems, which have shown both slight 
and substantial planar misalignmen t via observations o f 
disk orientation (jJensen et all 12004; IMonin et all 120061) . 
Probing the rotational evolution of intermediate separa- 
tion binaries (~1 - 10 AU) in the substellar regime will 
determine whether such trends hold at the lowest masses, 
with interesting implications for the formation and evo- 



lution of all types of binary stars. 

In this paper, we present projected rotational veloc- 
ity measurements for the components of a sample of 
tight, visual VLM binaries. The measurements of these 
spatially-resolved velocities are enabled by the W.M. 
Keck Observatory laser guide star adaptive optics (LGS 
AO) system, which provides high spa tial resolution ob- 
serva tions of optically faint targets ( Wizi nowich et all 
120061 ) . This study is the first to systematically examine 
the rotational velocities of individual VLM objects that 
reside in binary systems. In Section [2j we describe our 
observations and our method for extracting rotational 
velocities from high resolution spectra. In Section [3l we 
compare our measurements to those of single VLM ob- 
jects, and discuss the implications of our measurements 
for theories of binary star formation and evolution. We 
also discuss our results in the context of previously mea- 
sured radio emission from two of our sources. We sum- 
marize our findings in Section [4] 

2. DATA AND ANALYSIS 

2.1. Sample 

Our sample is comprised of eleven VLM binaries that 
were targeted as part of an ongoing program to measure 
their dynamical masses. These objects have been ob- 
served both astrometrically and spectroscopically since 
2006, and initial estimates of the ir orbital properties hay e 
been obtained from this dataset (jKonopackv et all2010f l. 
Their spectral types range from M7.5 to L7.5, and their 
separations range from 0'.'07 to 0'.'35. Because we are able 
to spatially resolve the components before obtaining high 
resolution spectroscopy (see section |2~2")) . our total sam- 
ple consists of 22 VLM objects. Tabic Q] summarizes the 
targets in our sample. 

2.2. Observations 

The eleven binaries were observed using the NIR spec- 
trograph NIRSPEC on Keck II 10 m (McLean et al 
2000) in conjunction with the facility LGS AO system 
(NIRSPAO). These observations, taken between 2006 
December and 2011 June, are described in detail in 
Konopacky et al. (2010). Briefly, we used the instru- 
ment in its high spectral resolution mode, selecting a slit 
0'.'041 in width and 2'.'26 in length in AO mode. We ob- 
served in the K band in order to obtain data in the CO 
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Table 2 

Log of NIRSPAO LGS K-band Observations' 1 



Target Name 


Date of 


A0V Star 


Exposure Time 


No. of 


Avg. SNR 


Avg. SNR 




Observation (UT) 


Standard 


(sec) 


Frames 


Primary 


Secondary 


2MASS J07464256+2000321AB 


2006 Dec 16 


HIP 41798 


1200 


4 


52 


44 




2007 Dec 04 


HIP 41798 


1200 


6 


72 


59 




2008 Dec 19 


HIP 41798 


1200 


6 


66 


56 




2009 Dec 09 


HIP 41798 


1200 


2 


83 


70 


2MASS J14263161+1557012AB 


2007 Jun 08 


HIP 73087 


1200 


1 


11 


33 




2008 Jun 01 


HIP 73087 


1200 


1 


50 


36 




2009 Jun 12 


HIP 73087 


1200 


4 


41 


29 


2MASS J17501291+4424043AB 


2008 May 31 


HIP 87045 


1200 


4 


48 


36 




2009 Jun 12 


HIP 87045 


1200 


6 


11 


31 




2010 Jun 07 


HIP 87045 


1200 


4 


36 


26 


2MASS J18470342+5522433AB 


2007 Jun 08 


HIP 93713 


1200 


4 


69 


60 




2008 Jun 01 


HIP 93713 


1200 


5 


69 


60 




2009 Jun 13 


HIP 93713 


1200 


3 


39 


36 




2010 Jun 07 


HIP 93713 


1200 


4 


53 


47 


2MASS J21402931+1625183AB 


2007 Jun 09 


HIP 108060 


1200 


4 


13 


28 




2008 May 31 


HIP 108060 


1800 


3 


58 


40 




2009 Jun 13 


HIP 108060 


1800 


2 


38 


26 




2009 Dec 10 


HIP 116611 


1800 


3 


56 


45 


2MASS J22062280-2047058AB 


2007 Jun 09 


HIP 116750 


1200 


3 


17 


39 




2008 Jun 01 


HIP 109689 


1200 


4 


54 


48 




2009 Jun 12 


HIP 109689 


1200 


4 


47 


44 




2010 Jun 07 


HIP 109689 


1200 


4 


29 


29 


GJ 569Bab 


2007 Jun 09 


HIP 73087 


900 


2 


89 


82 




2009 Jun 13 


HIP 73087 


900 


4 


86 


67 




2010 Jun 06 


HIP 73087 


900 


4 


111 


97 


HD 130948BC 


2007 Jun 09 


HIP 73087 


1200 


4 


13 


37 




2010 Jun 07 


HIP 73087 


1800 


4 


59 


52 




2011 Jun 18 


HIP 73087 


1800 


4 


52 


48 


LHS 2397aAB 


2007 Dec 04 


HIP 58188 


1800 


2 


68 


27 




2008 May 31 


HIP 61318 


1800 


3 


114 


41 




2008 Dec 19 


HIP 58188 


1800 


3 


83 


31 




2009 Jun 12 


HIP 61318 


1800 


2 


103 


33 




2009 Dec 09 


HIP 58188 


1800 


3 


105 


29 




2010 Jun 07 


HIP 61318 


1800 


2 


77 


31 


LP 349-25AB 


2006 Dec 16 


HIP 5132 


600 


4 


58 


45 




2007 Dec 04 


HIP 5132 


900 


1 


63 


58 




2008 Dec 19 


HIP 5132 


1200 


4 


105 


81 




2009 Jun 12 


HIP 5132 


1200 


4 


114 


98 




2009 Dec 09 


HIP 5132 


1200 


4 


121 


107 


LP 415-20AB 


2008 Dec 19 


HIP 24555 


1200 


4 


12 


32 




2009 Dec 09 


HIP 22845 


1800 


2 


62 


49 



a All data taken before December 2009 represents the same NIRSPAO-LGS data set presented in lKonopackv et all 1 12013) 



bandhead region (2.291 - 2.325 /im, order 33). Due to 
the dense population of lines in this region, our analysis 
for this work was done only in order 33, although the 
cross-dispersed data ranged from 2.044 - 2.382 ^m. 

The camera was rotated such that both components 
of each binary fell simultaneously on the high resolution 
slit, which is at an angle of 105.9° with respect to ver- 
tical. Typical observations consisted of four spectra of 
both components, each with 900-1800 second integration 
times, taken in an ABBA dither pattern along the length 
of the slit. On average, we achieved Strehl ratios between 
10-40% at K band, resulting in PSF c ore full width half 
maxim a of ~0'.'05-0'.'08. As discussed in lKonopackv et all 
(2010), this performance allowed us, in general, to obtain 
resolved spectra for binaries separated by at least ~0'/06. 

Table [2] gives the log of our spectroscopic observations, 
listing the targets observed, the date of observation, the 
number of spectra, the integration time for each spec- 
trum, and the average SNR achieved. Because the spec- 
tra are dominated by molecular features (see Figure [2]), 
we have chosen to estimate effective SNR by calculat- 
ing the average number of electrons per pixel in an ex- 
tracted spectrum and then assuming Poissonian noise on 



that average. We have verified that these estimates are 
roughly correct by calculating the SNR on small regions 
of our most rapidly rotating source, where all features 
are fairly smoothed out, and by using the properties of 
the NIRSPEC detector under the assumption that we 
are background limited at K baneful. Each target obser- 
vation was accompanied by the observation of a nearby 
A0V star to measure the telluric absorption. 

2.3. Data Reduction 

As described in iKonopacky et "aLl (|2010f ) . the basic re- 
duction of the NIRSPAO spectra was performed with 
REDSPEC, a software package designed for NIRSPEC0- 
Object frames are reduced by subtracting opposing nods 
to remove sky and dark backgrounds, dividing by a flat 
field, and correcting for bad pixels. Order 33 was spa- 
tially rectified by fitting the trace of each nod of AO cali- 
brators with third order polynomials, and then applying 
the results of those fits across the image. As these sys- 
tems are fairly tight binaries, cross-contamination can 

12 http:/ /www2. keck. hawaii.edu/inst/nirspec/sens. html 

13 http:/ /www2. keck. hawaii.edu/inst/nirspec/redspec/index. html 
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2.295 2.300 2.305 2.310 2.315 2.320 

Wavelength (urn) 

Figure 1. Example fits for two objects in our sample, GJ 569Ba 
and LP 349-25A. The spectra have been normalized and sepa- 
rated on the y-axis by a constant vertical shift for visual clarity. 
The black solid lines show the actual NIRSPAO data (order 33), 
which is not corrected for telluric absorption. The red dashed lines 
show the fitted template, which is a combination of a theoretical 
PHOENIX spectrum and a telluric absorption model. The i;sini 
values measured from these particular fits are 18 km s- 1 for GJ 
569Ba and 56 km s" 1 for LP 349-25A. For additional examples of 
fits fo r all sources in our sample, see Figure 4 of Konopackv ct al. 
pulOh . 

be an issue when extracting the spectra. We use a Gaus- 
sian extraction method to extract the spectra, fitting the 
trace with a variable FWHM in cases where the sep- 
aration was greater than 7 pixels and a fixed FWHM 
if less than 7 pixels. We subtract the results of this 
fit for one component before extracting the other (see 
iKonopackv et al1l2010l for more details). 

2.4. Determination of vsini 

The data set pres e nted here is the same set presented 
in IKonopackv et all (j2010T ) , except we now include three 
new epochs o f data from December 20 09, June 2010, and 
June 2011. In IKonopackv et ail (|2010ft we were primarily 
interested in the radial velocities, and hence orbital solu- 
tions, that could be derived from these spectra. Here we 
reanalyze these data having implemented two changes 
to our analysis in order to properly determine v sin i 
(Bai lev et al][2012T) . First, we now use the telluric lines 
present in the spectra of the A0V calibrator stars to mea- 
sure the instrumental profile for NIRSPAO. Secondly, we 
are now performing the convolution of theoretical tem- 
plates with a Gaussian kernel after putting all spectra 
on a log-linear scale. This makes the resolution constant 
across the entire spectral range, providing a more accu- 
rate measure of v smi. We describe our analysis in more 
detail below. Note that this reanalysis does not substan- 
tially impact our radial velocity estimates, which will be 
presented in a future paper. 

It has been demonstrated that the CO bandhead line 
depths and widths are primarily a function of temper- 
ature and the projected rotational velocity (t?sini) for 
VLM objects, respectively, with an additional mod erate 
dependence on surface gravity (Blake et al. 2007]) and 
metallicity. With some knowledge of the temperature 
of a given object and an allowance for unknown surface 
gravity and metallicity, v sin i measurements can be esti- 
mated from our extracted spectra. 

Our extracted spectra are not corrected for tel- 
luric absorption because these features provide a sta- 



ble reference for absolute wavelength calibration. Us- 
ing features that are naturally present in all spec- 
tra also allows us to accurately calibrate the instru- 
mental profile without the need to observe additional 
template sources. We therefore model each spectrum 
as a combination of a KP NO/FTS telluric spectrum 
(|Livingston fc Wallace! fl99ll ) and a synthetically gener- 
ated spe ctrum derived from the PHOENIX atmosphere 
models (jHauschildt et al.1 [199 9) . The model spectrum 
is parameterized to account for the wavelength solution, 
continuum normalization, instrumental profile (assumed 
to be Gaussian), vsini, and radial velocity. As men- 
tioned above, the instrumental profile is determined us- 
ing our A0V calibrator stars, which by design are a clean 
measure of the actual telluric spectrum. We hold the 
instrumental profile fixed while fitting our actual target 
spectra. The average resolution of our NIRSPAO data 
is ~10 km s _1 . The best-fit model is determined by 
minimizing the variance-weighted reduced \ 2 of the dif- 
ference between the model and the extracted spectrum, 
once this difference has been Fourier filtered to remove 
the fringing prese nt in NIRSPEC K-band spectra (see 
IBailev et al.l 120121 for more details) . This model there- 
fore provides our v sin i estimates. 

Each PHOENIX template is generated at a fixed tem- 
perature and surface gravity. The use of theoretical 
rather than observed templates has the advantage of in- 
troducing less template mismatch biases. Our main tem- 
plates for each source ha ve a temperature as measured in 
IKonopackv et al.l ((2010), a log(g) of 5.0, and solar metal- 
licity. Figure [1] shows example fits for two objects in our 
sample, one moderate rotator (GJ 569Ba) and one rapid 
rotator (LP 349-25 A) that have different temperatures. 
The figure demonstrates the location and morphology of 
the telluric features that are used for estimating the in- 
strumental profile and w avelength solution. Figure 4 of 
IKonopackv et all (|2010f) shows example fits for all ob- 
jects in our sample - we refer the reader to this work for 
further visual evaluation of our fitting technique. 

Statistical uncertainties are assigned by fitting each 
individual spectrum separately and taking the RMS of 
the values derived for each case. We also need to ac- 
count for systematic uncertainties due to both the tem- 
perature and surface gravity dependence of our spec- 
tra. We fit each spectrum with templates spanning 
±300 K in temperature and ranging from 3.0 - 4.5 dex 
in log(g). We also explored varying metallicity, using 
templates between ±0.25 dex of solar to fit our spectra 
based on metallicity measureme nts of low mass objects 
in the solar neighborhood (e.g.. IJohnson fc Appsl 120091 : 
iSchlaufman fc Laughlir] 120101 : iRoias-Avala et al.l 12010ft . 
We then use the spread in these values around our best- 
fit value as our systematic uncertainty, and add these in 
quadrature with our statistical uncertainties. We find 
on average that log(g) uncertainties add a 3 km s _1 un- 
certainty to the vsini, while temperature and metallic- 
ity uncertainties contribute an additional 1 km s _1 each, 
with lower values of temperature and metallicity yielding 
lower v sin i. The v sin i measured for each source at each 
epoch, along with the weighted average of all epochs, is 
given in Table |3] (all sources were observed at least two 
times). 

In order to confirm that our method returns the correct 
v sin i values, we obtained NIRSPAO observations of two 
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Figure 2. Single epoch examples of spectra for all objects in our sample (NIRSPAO order 33). All spectra have been normalized to the 
continuum and shifted to the same radial velocity. The spectra have been separated on the y-axis by a constant vertical shift, such that the 
absolute value of the flux is arbitrary. Sources are arranged in order of increasing rotational velocity and color-coded such that components 
in the same system have matching colors. Increasing rotational velocity dramatically impacts the morphology of the CO bandhead, which 
falls near the end of the K band at ~2.3/im. 



previously-measured M type stars. These objects, GL 
1245A (M5.5V) and G188-38 ( M4V), where targeted by 
several studies in the optical. iMohantv fc Basril (|2003[ ) 
measured projected rotational velocities of 22.5 ± 3.7 km 
s" 1 for GL 1245A and 29.4 ± 6.2 km s" 1 for G188-38. 
Other measurements for G1 88-38 include 36.5 ± 0.3 km 
s" 1 bv iDonati et alJ (j2006l) and 29.4 ± 1.4 km s" 1 by 
iDelfosse et al.l ( 19981 ). We performed an identical anal- 
ysis to that of our brown dwarf sample on these two 
mid-M stars, only using a higher temperature template 
that is more appropriate for these objects. We derived 
v sin is of 19 ± 3 km s" 1 for GL 1245A and 34 ± 3 km 



s _1 for G188-38, consistent with all the values from the 
literature. 

As an additional test, we can use "slow rotators" in 
our sample as templates, artificially spinning them up to 
estimate the vshii of other objects in our sample. Al- 
though this method suffers from template mismatch that 
is remedied by the use of theoretical atmospheres, it of- 
fers a further confirmation of our technique. We perform 
this analysis on one epoch of data for each binary, taken 
either in 2007 or 2008. We use the source with the lowest 
measured vsmi, 2MASS1847+55A, as our template for 
all other objects. This template, corrected for telluric ab- 
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Table 3 

v sin i Measurements (km s" 1 ) 



_L £11 J^U L 




Adopted 


2006 


2007 


2007 


2008 


2008 


2009 


2009 


2010 


201 1 


■ 

WClgiltCQ 




lypc 


T it (K} 3, 
1 nff (,J\I 


D c c c nib c r 


June 


D c c c nib c r 


]VI ay /June 


December 


June 


D c c c nib c r 


June 


June 


Average 


2MASS 0746+20A 


LO 


2205 


19 ± 5 




18 ± 5 




18 ± 5 




20 ± 3 






19 ± 2 


2MASS 0746+20B 


LI. 5 


2060 


33 ± 6 





32 ± 6 





32 ± 6 





34 ± 6 








33 ± 3 


2MASS 1426+16A 


M8.5 


2400 




6 ± 4 




6 ± 5 




7 ± 4 










6 ± 3 


2MASS 1426+16B 


LI 


2240 




8 ± 7 




12± 4 




11 ± 4 








11 ± 3 


2MASS 1750+44A 


M7.5 


2200 








9 ± 3 




10 ± 4 




8 ± 4 




9 ± 2 


2MASS 1750+44B 


LO 


2020 











10 ± 4 





10 ± 5 





12 ± 4 





11 ± 3 


2MASS 1847+55A 


M7 


2400 





3 ± 5 





4 ± 4 





9 ± 3 





8 ± 4 





7 ± 2 


2MASS 1847+55B 


M7.5 


2100 




5 ± 4 




5 ± 5 




9 ± 4 




9 ± 4 




7 ± 2 


2MASS 2140+16A 


M8.5 


2300 




13 ± 5 




12 ± 4 




11 ± 3 


16 ± 4 






13 ± 2 


2MASS 2140+16B 


L2 


2075 


— 


42 ± 7 


— 


37 ± 6 


— 


34 ± 6 


38 ± 7 


— 


— 


37 ± 3 


2MASS 2206-20A 


M8 


2350 




20 ± 4 




18 ± 3 




19 ± 4 




16 ± 6 




19 ± 2 


2MASS 2206-20B 


M8 


2250 




22 ± 4 




20 ± 4 




21 ± 4 




18 ± 6 




21 ± 2 


GJ 569Ba 


M8.5 


2000 




19 ± 3 








18 ± 3 




19 ± 3 




19 ± 2 


GJ 569Bb 


M9 


2000 




5 ± 5 








6 ± 5 




7 ± 4 




6 ± 3 


HD 130948B 


L4 


1840 




63 ± 8 












61 ± 7 


62 ± 6 


62 ± 4 


HD 130948C 


L4 


1790 




86 ± 8 












87 ± 11 


84 ± 12 


86 ± 6 


LHS 2397aA 


M8 


2180 






14 ± 3 


14 ± 3 


15 ± 3 


15 ± 3 


15 ± 3 


15 ± 3 




15 ± 1 


LHS 2397aB 


L7.5 


1350 






10 ± 6 


10 ± 6 


14 ± 6 


10 ± 8 


10 ± 6 


11 ± 8 




11 ± 3 


LP 349-25A 


M8 


2200 


56 ± 5 




50 ± 11 




59 ± 7 


55 ± 4 


54 ± 4 






55 ± 2 


LP 349-25B 


M9 


2050 


87 ± 6 




79 ± 11 




76 ± 10 


81 ± 4 


85 ± 5 






83 ± 3 


LP 415-20A 


M7 


2300 










40 ± 6 




41 ± 7 






40 ± 5 


LP 415-20B 


M9.5 


2000 










36 ± 4 




40 ± 7 






37 ± 4 



From Konopacky ct al. (2010). To be conservative, our analysis assumed a ±300 K temperature uncertainty for all objects 

3. DISCUSSION 



sorption using our observed AOV standards, was "spun 
up" to produce an artificial grid of spectra with v sin i 
between 5 and 100 km s . The grid was then cross- 
correlated with each object's spectrum, also corrected for 
telluric absorption, and we determined the value of v sin i 
that provided the maximum correlation. In all cases, the 
best value for v sin i found with this technique is within 
the uncertainties of the values given in Table [31 

Using these two independent checks, we are confident 
that our methodology is sound and that we are incorpo- 
rating the necessary uncertainties via our usage of mul- 
tiple temperature and log(g) templates. We do caution, 
however, that objects with particularly high values of 
■usini (>40 km s , four objects in our sample), though 
undoubtedly rapid rotators, might be subject to addi- 
tional systematic uncertainties not fully accounted for 
in our analysis due to greater sensitivity to properties 
associated with the instrument and technique, and this 
uncertainty may not be captured in the averaged values 
in Tabled 

We also estimate the lowest measureable value of v sin i 
in our spectra. To do this, we took our PHOENIX tem- 
plates and broadened them first to the correct instrumen- 
tal profile and then to different values of vs'mi. We also 
injected random Gaussian noise such that the templates 
would have SNR ~ 55, which is the average effective 
SNR per pixel for our data. We fit these spectra using 
the method described above. We find that the limiting 
value for which we could accurately measure v sin i is 3 
km s . Note that we are able to measure vsini be- 
low the instrinsic resolution of our NIRSPAO data due 
to our accurate measurement of the instrumental profile 
and theoretical templates that closely match our actual 
spectra. 

Figure [2] shows example spectra for all sources in our 
sample. These spectra, which we have corrected for tel- 
luric absorption for plotting purposes, are arranged in 
order of increasing wsini, demonstrating the effect of ro- 
tational velocity on CO bandhead morphology. 



This study represents the first measurement of compo- 
nent rotational velocities for a large sample of VLM bina- 
ries. The values presented in Table |3] show that ~80% of 
our sample are rapid rotators (wsini > 10 km s _1 ), and 
two sources, LP 349-25B and HD130948C, are among the 
fastest rotating VLM objects ever observed. In this sec- 
tion, we discuss the implications of these measurements. 

3.1. Comparison to v sini Measurements in the 
Literature 

In Figure |31 we plot our measured v sin is as a function 
of spectral type. We also include usini measurements 
from the literature, derived from seeing-limited obser- 
vations (IMohantv fe Basril 120031: IZapatero Osorio et al.l 
llMlRemers fe Basril2b^l20lbtlBlake et al.ll2010D with 
comparable sample, spectral resolution, and spectral 
type coverage to our observations. Our measurements 
are consistent with previous observations, which also 
find that VLM objects tend to be rapid rotators. In 
addition, our results are consistent with the trend of 
increasing v sin i with later spectral type. In a num- 
ber of cases, the values we measure are higher than 
those objects presented in the literature of a given spec- 
tral type. This is lik ely attributable to the mixed ages 
probed in this study. IReiners fe Basril (|2008) show that 
rotational velocity is a function of the age of the sys- 
tem, correlated with a spindown timescale that is driven 
by magnetic braking. In our sample, for instance, we 
find that the vs'mi measurements for HD 130948BC are 
higher than all previous measurements for mid-L dwarfs, 
and it has been proposed that this system is younger 
than the majority of the field population (~400 - 80 
Myr, IDupuv et all I2009at IMullan fe MacDonald|[20Toh . 
Indeed, the extremely rapid rotation of HD 130948BC 
m ay imply that the younger age of the system preferred 
by IMullan fc MacDonaldl lfeOlO) is more plausible. Our 
other extremely rapidly rotating system, LP349-25AB, 
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Figure 3. Projected rotational velocity (vsini) versus spectral 
type for each of the binary components in our sample (red cir- 
cles). Also plotted (open squares, uncertainties omitted for clar- 
ity) are previously measured values (seeing limited observations, 
so binaries ar e not resolved) of vs'mi for field VLM objects in 
the literature (Mohantv & Basri 2003; Zapatcro Osorio ct al. 2006; 
IReiners fc Basril200ai2010l:IBlake et al.120101 ). The values we mea- 
sure for our sources are consistent with other VLM objects, with 
our sources tending towards rapid rotation (vsini >10 km s —1 ). 

has also been prop osed to be quite young (~140 Myr, 
IDupuv et all [2010). On the whole, however, the mea- 
surements in our sample are fully consistent with the 
bulk population of objects previously observed. 

We also note that with these very rapid rotation rates, 
both LP 349-25B and HD 130948C are rotating at ~30% 
of their break-up speed. These large rotational velocities 
should also cause a high degree of rot a tional flattening. 
Using equation 6 of iBarnes fc Fortrievl (|2003|) . we derive 
that the ratio of the polar radii to the equatorial radii 
should be ~0.87. Given this level of oblateness, these 
objects might b e expected to exhibit me asureable lin- 
ear polarization (Seng upta fc Marlevll2010T) . In addition, 
the cooler temperatures at the equator due to gravity 
darkening may affect the spectral type measurements, as 
demonstrated by surfa ce imaging of rapidly rotating in- 
termediate mass stars (|Monnier et~aT1[2007l ). The rapid 
rotation also leads to a higher level of uncertainty in our 
radial velocity measurements for these extremely rapidly 
rotating sources. As is apparent in Figure [21 the major- 
ity of the CO bandhead features are basically smoothed 
by the rapid rotation, making anchoring these objects 
precisely in wavelength space quite challenging and giv- 
ing rise to the relatively hi g h rad ial velocity uncertainties 
given in IKonopackv et all ([2010D . 

Several sources in our sample were previously targeted 
in studies that did not resolve the components, but 
did measure vsini. LP 349-25AB, LH S 2397aAB, and 
2MAS S 2206-20AB were observed by IReiners fc Basril 
(2010). They obtained values consistent with ours for 
2MASS 2206-20AB and the primary component of LP 
349-25AB, but a very different value for LHS 2397aAB. 
We are not certain why our v sin i measurement of LHS 
2397aAB is different from these authors, but speculate 
that perhaps it can be attributed to the broadening of the 
lines due to the binary orbit in the unresolved spectra. 



20 40 60 80 100 

Primary vsin i (km s 1 ) 

Figure 4. The v sin i of each secondary component plotted against 
its primary. Sources with consistent velocities should fall on the 
dotted line. Five of our eleven systems show components with 
vsini that differ by >3<r. Of those five systems, four have sec- 
ondary components with higher velocities than their primaries. 



Uones et al.l (|199 6) measured the unresolved v sin i for LP 
415-20AB, also finding a result consi stent with our mea - 
surements (for both components). Bl ake et al.l (|2010l ) 
measured an unresolved vsini for 2MASS 0746+20AB. 
They obtain a nearly identical value to what we measure 
for 2MASS 0746+20B, although the combined light of 
the system should be dominated by the primary. How- 
ever, the flux ratio of ~1.4 at K band means the domi- 
nance is not extreme. 

3.2. Component v sin i Comparison 

We can also compare the rotational velocities of the 
components in each system to each other. The results of 
this comparison are shown in Figure [4j It is immediately 
apparent that a number of components in the same sys- 
tem have vastly different v sin i. There are five binaries in 
our sample that exhibit statistically-significantly (>3tr) 
differing v sini. These five systems, 2MASS 0746+20AB, 
2MASS 2104+16AB, GJ 569Bab, HD 130948BC, and LP 
349-25AB, have differences >10 km s -1 , or >30%. For 
all but GJ 569Bab, the secondary appears to be rotating 
more rapidly than the primary. For the other 6 systems 
in our sample, the velocities are consistent to within 2a. 
We also note that the consistency of these 6 systems, 
plus the generally correlated rapid or slow rotation of 
the 5 systems with different velocities, implies that the 
components of VLM binaries are not randomly paired in 
v sini. 

ISimon et al.l ()2006l ) noted that the lines of their 
spatially-resolved spectra were broader for GJ 569Ba 
than for GJ 569Bb. They measured v sin is of 25 km 
s -1 and 10 km s _1 for Ba and Bb, respectively, close to 
the values we derive here. The broadening of GJ 569Ba 
was postulated to perhaps be due to an unresolved third 
component rather than an intrinsic difference from GJ 
569Bb. Given that we see four other systems with differ- 
ing component vsini, it is even more plausible that this 
system d oes not have an unresolyed th ird component. In 
addition, Zap atero Osorio et al.l (|2004D measure v sin is 
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Figure 5. The component iisini (left) and absolute value of the difference between component vs'mi (right) as a function of the binary 
semi-major axis (a). All systems with components having v sini > 30 km s~ 1 and with significant iisini differences (>3<r) have a < 3.5 
AU. There also appears to be a rough trend of decreasing vs'mi and |A?;sini| with increasing a. The significance of these trends given our 
current data and relatively small sample is 2.3 - 2.9a for tisini vs. a and 1.4 - 2.4cr for |Ai>sini| vs. a. 



of 37 km s~ x and 30 km s _1 for Ba and Bb, thus also 
noticing a difference in the component values. The higher 
v sin i values potentially stems from their use of KI fea- 
tures in the J band that are known to be gravity sensitive. 
Although in many cases the rad ial velocity uncertaintie s 
for these sources are quite high (|Konopacky et al.ll2010r i. 
we do not see strong evidence for additional radial ve- 
locity variability in any of these discrepant systems that 
would point obviously to an additional unresolved com- 
ponent. 

We note that in the case of HD 130948BC, 
iMullan fc MacDonaldl (|2010t ) postulated that the com- 
ponents may not be rotating at the same rate, which we 
have now shown may be the case, although in contrast to 
their predictions, the secondary is likely the more rapid 
rotator and hence potentially the more magnetically ac- 
tive component. However, as we do not in fact know 
the true equatorial velocity (~v eq ) of either component, 
we cannot make any definitive statements about their 
models, although the rapid rotation of both components 
suggests that magnetic activity could be significant for 
both. See Section 13.41 for additional discussion of this 
system. 

The targets in this sample all ha ve orbital param- 
eter e stimates from previou s work s (iKonopacky et all 
20101 iZapatero Osorio eFall [20041: iSimon et all 120061: 
Dupuv et al.H2009allbH d. 120101) . We therefore explore the 



potential impact of the binary orbital properties on the 
consistency of usini. For the p urposes of this paper, w e 
use orbital elements derived in iKonopacky et all (|2010Tj . 
We looked for trends in both vsmi and Avsini as a 
function of all orbital parameters, and the only vari- 
able that produces a noticeable trend is semimajor axis 
(a). All components with vsmi > 30 km s" 1 and all 
five systems with significantly different v sin i have a < 
3.5 AU. In Figure we plot both the component vsmi 
and | Aw sin i | versus a. To assess the significance of the 
apparent trend of decreasing velocity and velocity dif- 
ference with increasing a, we use the Spearman's rank 



correlation coefficient. After accounting for the uncer- 
tainties in v sin i and a via Monte Carlo simulation, we 
determine using this metric that the significance of the 
trend in vsmi is between 2.3 and 2.9a and the signifi- 
cance in Avsini is between 1.4 and 2.4er. Because of 
the relatively small size of our sample, we are unable to 
explore the significance of this trend in greater detail. 
However, we note it here as a possible relationship of in- 
terest in the context of the discussion belo w. In addition, 
this tr end is similar to what was seen in Patience et al. 
(2002) , who observed that tighter binaries in a Per were 
rotating more rapidly than w ider binaries. Howeve r, such 
a trend was not observed by iBouvier et al.l (|1997t ) in the 
Pleiades. 

Given these result for the five "discrepant" systems, 
we are confronted with two possible scenarios. Either 
the binary components in some systems truly rotating 
at significantly different rates, or their rotation axes are 
inclined with respect to each other and possibly their 
orbital plane (or some combination of these two). We 
explore these possibilities below in sections 13.31 and 13.41 

3.3. Intrinsic Rotational Velocity Differences? 

It is possible that the binaries in our sample with "dis- 
crepant" vsmi have parallel rotation axes but differing 
rotational velocities due to intrinsic processes at work 
during either their formation or early evolution. Binary 
systems are generally thought to form via fragmenta- 
tion of a molecular cloud core or large circumstellar disks 
wherein small seeds are formed that eventually accrete 
more ma terial, form a disk, and achieve dynam ical stabil- 
ity (e.g.. lBonnell et al.lfl99lt iShu et al.lll990l ). Although 
the formation mechanism for VLM objects is still an open 
question, a fragmentation origin is certainly plausible for 
these objects. Simulations of core fragmentation origi- 
nally assumed that the rotation axes of the binary seeds 
were aligned with the rotation axis of the core (|Batd 
1 1997T ) - In these simulations, it was shown that the prop- 
erties of the specific angular momentum of the accreting 
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Figure 6. The measured and expected Avsini (secondary mi- 
nus primary) for each binary based on the spectral type of the 
components. The expected value is derived using the relationship 
presented in Rcincrs & Basri (2008) for vs'mi versus spectral type 
at an age of 2 Gyr. We assume an uncertainty in derived usini 
for each source of 20%, which gives uncertainties in Ausini of ~7 
to 13 km s . While this relationship can account for the Avsini 
of 2MASS 0746+20AB and 2MASS 2140+16AB (marginally), the 
values of GJ 569Bab, HD 130948BC, and LP 349-25AB are all 
>1.6ct discrepant. 

material onto the protobinary had a substantial impact 
on the properties of the binary. For a binary to grow to 
a mass ratio of about 1, as in the case of most of the 
objects in our sample, requires a higher specific angu- 
lar momentum for the accreting material, which in turn 
leads to a higher accreti on rate for the secondary compo- 
nent than the primary (|Batdll997D . The fraction of this 
angular momentum that is converted into orbital angu- 
lar momentum versus spin angular momentum depends 
on the size of the accretion radius. It has been shown, 
however, that in cases where significant spin increase is 
achieved from accretion, the primary component tends 
to increase more than the secondary component (jBatel 
U997tlArtvmowiczlll983l) . which is not in agreement with 
most of our "discrepant" systems. In addition, if these 
objects formed circumstellar disks after the initial accre- 
tion phase, they may be subjected to braking through 
magnetic coupling with the disks, as discussed below. 

As has been demonstrated in other works and dis- 
cussed earlier, there is a clear evolution of rotational ve- 
locity as a function of age for VLM o bjects, much like 
what is observed for higher mass stars (|Reiners fc Basril 
2008). It has been suggested that during the pre- 
main sequence phase of evolution, rotation is likely reg- 
ulate d by magnetic coup l ing to a circumstellar disk 
(e.g., lEdwards etahl 119931: iKundurthv et alj I2006Q . al- 
though this mechanism is still a matter of debate (e.g., 
IStassun et all 119991: iNguven et al.l 120091). O ncc the disk 
dissipates in ~1-10 Myr (| Strom et alj 11993ft . a star or 
VLM object will speed up due to its continued contrac- 
tion. It has been shown that circumste llar disks exist 
around components of b inary stars (e.g., iMcCabe et alj 
2006; Cic za et ail 12009). If component disks dissipate 
on different timescales, one object will begin to speed 
up sooner than the other. This could potentially lead to 



differing rotational velocities in spite of coevality. Obser- 
vational support for this possibility exists via measure- 
ments of disks aroun d T Tauri binaries. For instance, 
IMcCabe et"aT1 (|2006ft identify six T Tauri binary sys- 
tems that have component disks in different phases of 
evolution, with the secondary tending to have the more 
evolved disk (albeit for slightly higher mass objects than 
observed in this work). However, it is unclear if a dissipa- 
tion timescale difference, which would likely be at most 
around 10 Myr, is alone sufficient to generate a large 
rotational velocity difference in these systems. Further- 
more, the tight physical separations of the binaries in our 
sample should have impacted the formation and survival 
of any circumstellar disks, leading to both components 
having disks that last for at most 1 Myr (e.g.. lCieza et al.1 
I2009t IDuchenel 12010ft . The truncated disk survival time 
in close binaries could be related to the tentative trend 
we see in v sin i as a function of a in Figure [5j 

After the disk dissipates, VLMs likely follow a slightly 
modified form of the wind braking law that has been 
shown to reproduce the rotation versus age correlation 
seen amongst higher mass objects. The timescale for 
the braking of these objects is significantly longer than 
for high mass objects. In all cases, the rotational speed 
seems to be clearly a function of mass, which in a co- 
eval system is correlated with spectral type. It could 
be postulated, therefore, that the reason we observe 
some secondaries rotating more rapidly than primaries 
is due to thei r lower mass/later spectral type. Indeed, 
IReiners et al.l ()2007ft identify this as the cause of the 
v sin i differences in the components of the triple system 
LHS 1070. In order to determine whether this is a viable 
cause of the differences in v sin i we observe for our bina- 
ries, we attempt to roughly assess the expected vsini as 
a function of spectr al type for a given age. We use the 
relationship given in IReiners fc Basril (|2008ft , which uses 
the wind braking law to describe the rotational velocity 
evolution of VLM objects, for an age of 2 Gyr. We choose 
this age because, although the ages of the sources in our 
sample are generally unknown, it is roughly the average 
age of the systems in our sample that are thought to be 
relatively young (<500 Myr) and those that are thought 
to be the age of the field (possibly as old as ~5 Gyr). Us- 
ing this relationship, we calculate the expected vsini for 
each binary component based on its spectral type, and 
then compute the expected difference between the two. 
We also assume an intrinsic scatter at a given spectral 
type to account for some of the spread seen in measure- 
ments of rotation rates in clusters with known ages (i.e., 
iTerndrup efaLI 120001: llrwin et al.1 12009ft . We assume a 
conservative intrinsic scatter per object of 20%, which 
gives an uncertainty of between ~7 and 13 km s _1 in 
Avsini. Because of the fairly substantial intrinsic scat- 
ter assumed here, our choice o f usin g the 2 Gyr rela- 
tionship from IReiners fc Basril (|2008ft has little impact 
on this comparison because values of Av sin i from their 
2, 5, and 10 Gyr relationships are all consistent given 
these uncertainties except in the case of LHS 2397aAB. 
For this system, older ages predict an even greater dif- 
ference in Avsini, whereas we observe no statistically 
significant difference in v sin i in this system. The only 
assumption that could potentially yield Av sin i consis- 
tent with our measurement is a very youn g age, which 
is most likely not the case for this system (|Freed et al.1 
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Figure 7. Joint PDFs between the required relative inclination of 
the components of the five "discrepant" binaries and their assumed 
v e( j • v e9 is determined by assigning the more rapidly rotating com- 
ponent an inclination sampled from a distribution uniform in cos(i) 
(ignoring the known orbital inclination) and is assumed to be equal 
for both components. The white circle on each PDF denotes the 
peak probability. All the peaks fall between 40 and 60°. Small 
values of Ai generally require very high v e9 . The figures only ex- 
tend to 300 km s _1 , which is roughly the break up speed for these 
objects. 



I2003I: IDupuv et afl 12009b). We plot the values we de- 
rive for expected Av sin i in Figure [SJ along with our 
actual measurements for each binary. While this ef- 
fect, given our assumptions, could explain the differences 
in the component velocities in 2MASS 0746+20AB and 
2MASS 2140+16AB (marginally), it gives results incon- 
sistent with our measurements for GJ 569Bab (1.6<7 off), 
HD 130948BC (1.6a off), or LP 349-25AB (2.9ct off)! 
And, as mentioned, it predicts a large velocity difference 
for the components of LHS 2397aAB, which we do not 
observe, ft is also worth noting that GJ 569Bab, HD 
130948BC, and LP 349-25AB have all been postulated 
to be fairly young (<700 Mvr.lDupuv et alj|2009al [2011 
iSimon et alboOrj iZapatero Osorio et al.ll2006l ). 

3.4. Mutually Inclined Rotation Axes? 

If we instead assume that the binary components with 
similar masses must have similar or the same rotational 
velocities, then their rotation axes must be inclined with 
respect to one another in systems with differing usini. 
To determine the most likely value of the relative inclina- 
tion of the spin axes in these five systems, we ignore for 



a moment the known orbital inclination and perform a 
Monte Carlo simulation assuming one component has an 
inclination sampled from a distribution that is uniform 
in cos(z). We select a value of vs'mi for each compo- 
nent from a Gaussian distribution defined by our mea- 
surements and uncertainties, and assign the randomly 
sampled inclination to the faster rotating component to 
determine the equatorial velocity (v eq ) for the system. 
We then calculate the inclination required for the other 
component to have the same v eq . We perform this exer- 
cise 100000 times to derive probabilities, which we plot 
in Figure [7J Note that there are two possible values for 
the component with "unknown" inclination that will give 
v eq due to a 180° ambiguity (axis pointing up or down 
from our line of sight). Figure [7] shows joint probability 
density functions (PDFs) between v eq and the absolute 
value of the difference in the inclination of the two com- 
ponents. The white dot ovcrplottcd on each distribution 
denotes the location of peak probability. These peaks are 
found between ~40 and ~60 degrees in Ai and fall close 
to a v eq that equals the v sin i of the faster rotating com- 
ponent since a randomly oriented object is more likely 
to be observed edge-on than pole-on. The figure also 
demonstrates that the only configuration that can main- 
tain both the same inclination and v eq is one in which 
both components rotate close to break up speed (~300 
km s _1 ) and are observed almost pole on. This is an 
extremely unlikely configuration and can safely be ruled 
out as an explanation for all the systems with discrepant 
vsini. 

It is possible that some level of rotation axis misalign- 
ment is natural in these binary systems and could be 
represented by a Gaussian distribution centered on some 
"typical" misalignment. We can use our sample to as- 
sess this toy model assuming the components have the 
same v eq . In order to determine the value for the typi- 
cal relative inclination and its la spread, we performed 
another Monte Carlo simulation in which we sampled a 
relative inclination for each system from the distribution 
allowed by our vsini measurements and uncertainties. 
In this case, however, we assumed that one component 
was aligned with the orbital plane, sampling from a dis- 
tributio n allowed by the measu rements of orbital incli- 
nation ([Konopackv et all 120101 ) . These inclination dis- 
tributions tend to avoid cases with very high v eq and 
therefore are more realistic. As before, we assigned the 
more rapidly rotating component the sampled inclina- 
tion, and calculated the inclination required for the other 
component to have the same v eq . We then generated 
many (100000) distributions of relative inclinations for 
our sample, which we fit with a simple Gaussian model 
to find the peak and full width half maximum. We per- 
formed the simulation for two cases, one in which we 
chose the smaller of the two allowed relative inclinations 
(due to the 180° ambiguity), and the other chosing at 
random either the smaller or the larger of the the al- 
lowed relative inclinations. In the first case (small an- 
gles) we find that the preferred distribution has the form 
18 ± 26° and in the second case (any angle) the distri- 
bution has the form 25 ± 60°. The implication of this 
is that if we are truly probing a Gaussian distribution of 
relative inclinations with the expectation of a few many- 
sigma outliers, the average inclination and spread must 
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Figure 8. Top two panels: Eccentricity oscillations and re- 
gression of the ascending node (in the reference frame of the A-BC 
orbit) for HD 130948BC due to the presence of HD 130948A. Based 
on the obse rved and mo deled properties of the system, it could be 
undergoing Kozai (1962) oscillations on a ~ 10 yr period. Bot- 
tom three panels: Evolution of the angular difference in the spin 
vectors (Aig c), rotation axis inclination to the line of sight, and 
a possible resulting time history of v sin i of HD 130948B (solid 
line) and C (dashed line) due to the presence of HD 130948A. This 
effect could explain our v sin i measurements and do so within the 
expected lifetime of the system (~500 Myr.) 

be quite substantial in order to describe for our sample. 

The idea of forming objects in which the rotation axis is 
inclined with respect to the orbital axis has been explored 
in great detail recently due to the discovery of planets 
orbiting in a plan e that is misaligned wi t h the stellar 
rotat ion axis (e.g., iWinn et al.ll2009|. | 2010fc iTriaud et al.l 
I2010H . For instance. iBate et al.l (|2010f ) investigated the 
accretion history of stars forming in a turbulent cluster 
environment. Given variable accretion rates and mate- 
rial being accreted from different directions, it is possible 
to impact both the rate of rotation and the axis of ro- 
tation, leading to rotation axes th at are mis a ligned with 
circumstellar disks. Alternatively, lLai et all l)201ll ) have 
explored how the interaction between the magnetic field 
of a young star and its circumstellar disk can effectively 
push the stellar rotation axis out of alignment due to a 
magnetic warping torque. Such scenarios could poten- 
tially cause binary stars to have inclined rotation axes. 

A third possibility is s ecular torqu es due to a third stel- 
lar companion. Indeed, iHald (|1994f ) gave a general rule 
that binary components closer than 40 AU have agree- 
ment in spin orientation, but mentioned that rule is bro- 
ken for systems with a third body. In particular, consider 
a binary star that is born with each of its component spin 



vectors aligned with the orbital angular momentum. If 
the orbital orientation remains fixed throughout the bi- 
nary's lifetime, there would be no torque to cause spin 
precession, and the components would remain aligned 
with the orbit. However, if a companion star (777,3) is 
introduced on an external, non-coplanar orbit of period 
P ou t and inclination i relative to the binary (mi-7772) of 
period Pj n , it will cause orbital nodal precession in the 
binary on a secular period: 



P out ™i 



+ m 2 + 777-3 



3 cos i Pi, 



7773 



(1) 



(|Kiseleva et al.|[T998f ). In response to this misaligned or- 
bit, the spins will precess as well, on a period: 



process, 1 



777 2 fel \Rl 



; (1 



=2^3/2 



COSIp 



(2) 



(jEggleton fc Kiseleva-Eggletonl 12001) for star 1, where 
fc2.i is the apsidal motion constant /c 2 for star 1, C% is its 
normalized inertia (= I /MR 2 ), Pi is its spin period, ip 
is its spin obliquity, and a and e are the semi-major axis 
and the eccentricity of the orbit. An analogous period 
for m 2 need not be the same if the masses and radii are 
not quite equal. Thus we have the opportunity to begin a 
binary aligned, and due to slight differences in precession 
rate, open up a large inclination difference between the 
component spin axes. Thus, for the spins to become 
misaligned with one another, we require 



precess,! ^precess, 2 



1 



system age 



C-i) 



However, if the two spin precession periods are differ- 
ent but much shorter than the orbit precession period 
(equation [1]) , then the spins will just track the orbit ori- 
entation and will not become misaligned from either the 
orbit or each other. Therefore, we also require: 



Psec 777(73? (Xp reC ess,l , Pprccess,2) 



(4) 



In our sample, we have two binaries with similar spec- 
tral types (GJ 569Bab and HD 130948BC) that are part 
of known triple systems. We have therefore produced a 
numerical demonstration of this effect for HD 130948BC, 
which has component vsini values that differ, despite 
similar luminosities and spectral type. Although the 
mass ratio of this system has no t been direct l y mea- 
sured, we can follow the method of IDupuv et al.l (|2009al) 
and use a system age estimate and the component bolo- 
metric luminosities to estimate individual masses. This 
method predicts a slight mass difference due to the fac- 
tor of ~1.1 difference in the luminosities. While the 
components should have very similar radii now, this 
was not the case when they were younger if they do 
indeed have slightly different masses. Because the ra- 
dius is such an important factor in equation [2] we as- 
sume for the purposes of this demonstration the average 
radii for these objec ts over the a p proxim ate age of the 
system (~500 M yr. IDupuv et all i2009aD ): (Ri,R 2 ) = 
(0.158, 0.150)i?@ (|Chabrier e t al. 2000, DUSTY models), 
and corresponding mas s es (Mi,M 2 ) = (0.0 6, 0.05)Af Q 
(jKonopackv et all 120101: IDupuv et al"1l2009a|) . For the 
sake of this argument, we set Pi = 1.92 hours and 
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Pi = 1.82 hours (i.e., v eq = 100 km s _1 ), so that the 
observed difference in v sin i will be due to inclination 
differences rather than rotation rate differences. Last, 
we chose apsidal motion constants fc2,i = &2.2 = 0.175 
and normalized mom ent of inertia C\ = = 0.108 
(jLeconte et al.l 120 111 ). We integrated the binary orbital 
parameter s (starting with the obse rved a — 2.19 AU and 
e = 0.16, iKonopackv et all I2010T) in the gravitational 
potential of the third star (m 3 = 1.1M , a = 50 AU, 
e = 0.7 in an orientation i = 80°, uj in — 45°), using a sec- 
ular code that uses the quad rupole approximation for the 
dynamics of the three stars dFabrvckv fc Tremaindl2007l : 
Egglc ton fc Kiseleva-Eggletonl 120011^ In Figure |5] we 
show the results. On the short period T sec « 10 4 yr, the 
inner orbit processes and goes through eccentricity oscil- 
lations, i.e., iKozail (|1962f) cycles. On the longer period 
T precess w 4 x 10 s yr, the spins precess, and on an even 
longer period, 1/(1/T prcccss ,i - 1/T prcccss , 2 ) « 1 x 10 9 yr, 
the precession rates mix in phase (one gains by n, due 
to the slightly different masses and radii) . The system 
is likely as young as 500 - 70 Myr (|Dupuv et alj|2009at 
iMullan fc MacDonaldl l2010f ). Therefore it is plausible 
that the differential precession only recently opened a 
measurable angle between the rotation axes of these two 
components. In any case, this exercise suggests that no 
matter the true masses and radii of the components, 
Kozai cycles are quite possibly occurring in this sys- 
tem, and the differences in component v sin is are pos- 
sibly just one manifestation of the precession of the in- 
ner binary. The rapidity of the oscillation period implies 
that the precession of the inner orbit and/or its eccentric- 
ity oscillation may be measurable; e.g., the eccentricity 
currently has a precisio n of ~0.01 ()Dupuv et al.l [2009a; 
IKonopackv et all 120101 ) and this numerical integration 
shows eccentricity changes of 0.01 per decade. 

We have also performed an analogous simulation for 
the G J 569Bab system, which exhibits v sin i differences 
of only 13 km s _1 as opposed to the 24 km s _1 differ- 
ence in HD 130948BC. Using the orbital parameters of 
GJ 569Bab, which has a tighter separation (~1 AU), we 
find that the precession period for the system is much 
more rapid, on the order of 0.1 Myr. It is therefore sub- 
stantially easier to have precession occur for this sys- 
tem, making it plausible that secular perturbations are 
responsible for the v sin i difference in this system as well. 

Three of the five systems with d iffering v sin is are 
not part of known triple systems. lAllen et al.1 ()2007f ) 
performed deep imaging around 2MASS0746+20AB and 
2MASS2140+16AB and found no comoving companions 
between 40-1000 AU down to a mass limit of ~0.05 M . 
Although we are unaware of similar deep imaging for 
LP 349-25 AB, current all-sky surveys do not reveal any 
bright sources within 1000 AU. This does not, however, 
rule out the possibility that these sources were previously 
members of higher order multiple systems, or had an in- 
teraction with an unrelated object. We also note that to 
our knowledge, none of the sources in our sample with 
consistent usini have additional companions. 

3.5. Implications for Radio Observations 

Two of our targets (2MASS074 6+20AB and LP 349- 
25AB) are known rad io sources (|Antonova et al.l 120071 : 
iPhan-Bao et al.ll2007D . The components of these bina- 



ries exhibit rapid rotation (>19 km s _1 ). Both systems 
also have components with different wsini. This has in- 
teresting implications for determining which of the bi- 
nary components is emitting in the radio, as the radio 
measurements are typically unresolved. For instance, 
iBerger et all (pool ) used a v sin i for 2MASS 0746+20AB 
from unresolved measurements to derive a radius for the 
radio emitting component, which they assumed to be the 
primary. Using the spatially resolved measurements of 
v sin i for these objects rather than an unresolved value 
(which at 27 km s _1 is nearly exactly the average of the 
spatially resolved values of 19 and 33 km s _1 that we 
obtain here) gives slightly different values for the radius. 
We can now derive the predicted radius for each com- 
ponent without assuming which is the emitting source 
and determine which gives a m ore plausible r esult. If we 
make the same assumption as IBerger et al.l ([2009) that 
the rotation axis is pe rpendicular to the orbit al plane, 
which was updated by IKonopackv et al.l (|2010l ) to have 
an inclination of 138.2 ± 0.5°, we derive velocities of 29 
± 3 km s" 1 and 50 ± 3 km s _1 for the primary and sec- 
ondary, resp e ctivel y. Using the rotational period from 
IBerger et all (|2009D of 124.32 ± 0.11 minutes gives a ra- 
dius of 0.050 ± 0.005 R Q for the primary and 0.085 ± 
0.005 R Q for the secondary. Assuming the total system 
mass and preliminary estima tes of the component mass 
from IKonopackv et all (|2010l ) that the mass of the pri- 
mary is between 0.08 and 0.1 M while the secondary 
is between 0.05 and 0.07 M Q , evolutionary models pre- 
dict radii of 0.103 - 0.125 R© for for the primary and 
0.094 - 0.096 R© for the secondary. The implied discrep- 
ancy with the models is therefore 10cr for the primary 
and 1.8cr for the secondary. Given this, in conjunction 
with the unphysically small radius predicted for the pri- 
mary, we postulate that the secondary is the source of the 
radio emission, which may be related to its more rapid 
rotation. Although the radius of 0.085 ± 0.005 R Q is 
consistent to within the uncertainties of the value from 
IBerger et al.l ()2009l ). the resulting overprediction of the 
radius by evolutionary models becomes slightly less se- 
vere because if the emission is from the secondary, the 
models predict it will have a smaller radius. In either 
case, this analysis highlights the importance of obtain- 
ing fundamental parameters of binary components indi- 
vidually rather than from unresolved measurements. In 
addition, the results from our sample overall suggest that 
the assumption of rotational axis alignment with the or- 
bital plane should be treated with much greater caution. 

We can also speculate that the secondary component 
of LP 349-25AB is likely the radio source in the system, 
since it is rotating extremely rapidly. However, it is not 
improbable that both components are radio sources. Fu- 
ture VLBI observations will probe the true origin of the 
radio emission in both of these system (G. Hallinan et al., 
in prep) . Furthermore, our results provide a guideline by 
which sources previously unobserved in the radio should 
be targeted. Binaries that have at least one component 
with v sin i > 30 km s _1 would make excellent candidates 
for observation. 

4. SUMMARY 

Using the combination of high spatial and high spec- 
tral resolution afforded by the Keck II LGS AO system, 
we have measured component v sin is for a sample of 1 1 
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VLM binaries. Among the 22 objects measured, 80% 
are rapid rotators (vsini >10 km s -1 ), consistent with 
previous measurements for VLM objects. We found that 
5 of the binaries surveyed had components with usini's 
that differed by >3cr. We explored potential causes for 
these differences, which must stem either from intrin- 
sic velocity differences or from mutually inclined rota- 
tion axes, or a combination of both. Our analysis shows 
that perhaps both explanations are required to explain 
these five binaries, with two (2MASS 0746+20AB and 
2MASS 2140+16AB) having explainable intrinsic differ- 
ences due to differing spectral types, while the other three 
(GJ569Bab, HD 130948BC, and LP 349-25AB) likely re- 
quire a different explanation. We looked at the binary 
HD 130948BC as an example of a case in which secular 
torques are most likely causing the orbit and component 
spin vectors to evolve. This binary, in which Kozai oscil- 
lations are possibly at work, is an example of the impact 
of dynamical evolution on VLM binaries. One other bi- 
nary in our sample with differing component velocities, 
GJ 569Bab, is part of a triple system and a similar anal- 
ysis shows it too may be displaying the impact of sec- 
ular torques. For LP 349-25AB, which is not part of a 
known triple system and spectral type differences do not 
seem to account for its vastly different v sin is, the ex- 
planation may rest with past dynamical encounters or 
perturbations by higher mass objects. Our results also 
have implications for the previous measurements of radio 
activity in this system and 2MASS 0746+20 AB. We sug- 
gest that the secondary components are more likely the 
radio emitting sources in these systems. This stresses the 
importance of measuring fundamental parameters of bi- 
nary components individually rather than bootstrapping 
from unresolved measurements. 

Continued monitoring of these systems will improve 
the precision to which parameters such as component 
mass and radius are measured, allowing for the correla- 
tion of the properties with rotational velocity. This will 
provide a new handle on the ages of these objects. Fur- 
ther, the objects that we found to be rapid rotators that 
have not yet been surveyed for radio emission are ideal 
targets for observation. 
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